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ABSTRACT 



Context. Dark energy is now one of the most important and topical problems in cosmology. The first step to reveal its nature 
is to detect the evolution of dark energy or to prove beyond doubt that the cosmological constant is indeed constant. However, 
in the standard approach to cosmology, the Universe is described by the homogeneous and isotropic Friedmann models. 
Aims. We aim to show that in the perturbed universe (even if perturbations vanish if averaged over sufficiently large scales) the 
distance-redshift relation is not the same as in the unperturbed universe. This has a serious consequence when studying the 
nature of dark energy and, as shown here, can impair the analysis and studies of dark energy. 

Methods. The analysis is based on two methods: the linear lensing approximation and the non-linear Szekeres Swiss-Cheese 
model. The inhomogeneity scale is ~ 50 Mpc, and both models have the same density fluctuations along the line of sight. 
Results. The comparison between linear and non-linear methods shows that non-linear corrections are not negligible. When 
inhomogeneities are present the distance changes by several percent. To show how this change influences the measurements of 
dark energy, ten future observations with 2 % uncertainties are generated. It is shown the using the standard methods (i.e. under 
the assumption of homogeneity) the systematics due to inhomogeneities can distort our analysis, and may lead to a conclusion 
that dark energy evolves when in fact it is constant (or vice versa). 

Conclusions. Therefore, if future observations are analysed only within the homogeneous framework then the impact of inhomo- 
geneities (such as voids and superclusters) can be mistaken for evolving dark energy. Since the robust distinction between the 
evolution and non-evolution of dark energy is the flrst step to understanding the nature of dark energy a proper handling of 
inhomogeneities is essential. 

Key words, cosmology - dark energy - cosmological parameters - large-scale structure of Universe 



1. Introduction 

The most important problem in contemporary cosmology 
is the nature of dark energy. There are many ongoing 
attempts to explain the nature of dark energy. The 
simplest explanation assumes that dark energy does not 
evolve and is modeled by only one single constant — 
the cosmological constant. In this scenario, the cosmo- 
logical constant is considered as one of the fundamental 
constants of nature. However, most scientists are not 
satisfied with this type of solution and continue to inves- 
tigate the nature of dark energy. Most o f these attempts 
invoke new phys i cs: f(R) cosinologies (IBuchdahl 197C : 
Starobinskv 1980l: iKerner 1982t iBarrow fc Ottewill 1983 : 
Barrow fc Cotsakis 19881: iLi fc Barrow 2007 



Peebles fc Rat ra 2003': 'Copeland et al. 2006[) . phan 
tom fields (.Caldwell et al. 200 3) , holog ra phic co s 



mology (iKaloper fc Linde 19991 : 
Simpson 2007h . or accounting for 
such as Cha pl ygin gas 



Li 2004 iHsu 2004 : 



exotic substances 
dBertolami et al. 2004 : 



Bento et al. 20031: iDebnath et al. 20041 ). There are 
also alternatives which replace dark energy with an 
i nhomogeneity of a radius of app rox imately one Gpc 
(ICelerier 2000l: lleuchi et al. 2002t lAlnes et al. 200 ' 



Bplcjko 2008.; 



Alnes fc Ama rzguioui 2006t Enqvist fc Mattsson 200 



Garci'a-Bellido fc Haugb0llc 2008 : 



Kovama 20081 : iMoldenhauer et al. 20101) . brane 

word models (IRandall fc Sundrum 1999 ; 

Dvah et al. 200d 



Deffavet et al. 200 it lAlcaniz 2002 : 



Clifton et al. 20081: iBoleiko fc Wvithe 20091: 

Celerier et al. 2010 ) (for a review and explicit exam- 
ples see Bolejko et al. 2009) — these, exceptionally, do 
not require any new physics. In most of these alternatives, 
an effective dark energy equation of state evolves with 
time. Therefore, from the observational point of view, the 
most important discovery would be to detect any possible 



. . , ^^^^^^ iixi^v^i v,i^iiv, v^i^vjv^ . J .vv^.^iv* ^^^^^^ "-^V Jr^^'^'J^^^v^ 

Sahni fc Shtanov 20031: IMaia et al. 20051) . scalar fields changes in dark energy over the period of evolution of the 



(|Saini et al. 2000t ISahni fcWang 2000t IChiba et al. 2000t 
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Universe or to prove beyond doubt that the cosmological 
constant is indeed constant. 

The evolution of dark energy manifests itself in its 
equation of state. For the perfect fluid , the conservation 
equations T'^P.p = imply (jEUis 197l[ ) 



dpDE 

dt 



QipDE +PDe) = 0, 



where p is energy density, p is pressure, and Q is the ex- 
pansion parameter (in the Friedmann models Q — 3H; H 
is the Hubble parameter). Therefore, the only case when 
dark energy does not evolve is when poE = —POE- Thus, 
if only w = pde/pde differs from —1 then dark energy 
evolves. 

However, as long as dark energy is not directly detected 
any measurement of w is always indirect and relies on the 
assumption that some relations hold. In most cases, it is 
assumed that the distance-redshift relation is as in the 
Friedmann model. This paper shows that in the perturbed 
universe, even if {5p) = the distance relation is not the 
same as in the unperturbed universe. Therefore, the study 
of dark energy within the homogeneous framework may 
strongly bias the results of these analyses. 

Since the distance relation in its general form is a par- 
tial differential equation, to make the reader more aware 
of some features of differential equations a simple exam- 
ple with a harmonic oscillator is presented in Sec. [21 In 
Sec. 131 the relation for the distance in the general inho- 
mogeneous case is presented and an analogy with the har- 
monic oscillator solution is discussed. Sec. SI discusses the 
consequences of evolving dark energy. In Sec. [SI an exam- 
ple of the Swiss-Cheese model is presented to shown that 
small-scale inhomogeneities if not analysed properly can 
be mistaken for an evolving dark energy. 

2. Harmonic oscillator 

We now discuss some basic features of the harmonic oscil- 
lator. Although this material is most probably well know, 
it is instructive to review it here, as there is a important 
relation to cosmology, namely to the distance relation. The 
evolution of the harmonic oscillator is given by 



0. 



(1) 



X + UJqX 

The solution is a combination of sinujot and coswo^- We 
refer to this solution as the unperturbed solution. If some 
additional force is present (for simplicity we assume that 
the force is sinusoidal) then 



X + ujqX — Fq sincLit — 0. 



(2) 



The solution of this equation is a sum of the unperturbed 
solution and the steady-state solution given by 



Xs{t) 



Fn 



■ s'mujt. 



UJQ 



As seen, in this case, the driving force introduces oscil- 
lations of frequency uj around the unperturbed solution. 
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Fig. 1. The solutions of eqs. © - ©■ Sohd line: the un- 
perturbed solution, equation dH); dotted line: the solution 
of equation ^\ dashed line the solution of equation ([3]). 

However, if we depart from this highly idealistic scenario 
the solution changes. We assume that F is proportional 
to X (i^ = xFq sinwt), then 



X + ^(wq — -fo sinwt) = 0. 



(3) 



In this case, the solution does not oscillate around the 
unperturbed solution. This is presented in Fig. [J which 
presents the solutions of Eqs. (UJ - (O- The parameters 
used here are = 1, w = 15, and Fq = h. The unper- 
turbed solution, Eq. (jl}, is given by the solid line. The 
dotted line presents the solution of Eq. ([2]) - as seen, it 
oscillates around the unperturbed solution. The solution 
of Eq. ([3]) is given by the dashed line: it does not oscil- 
late around the unperturbed solutions, even though the 
average of F over the period vanishes, i.e. {F)t = 0. 

The above conclusion is rather obvious, although we 
show there is a close relation between the above exam- 
ple and the cosmology. We demonstrate that the distance 
in the perturbed universe (even if = 0) does not os- 
cillate around the homogeneous value. This has a strong 
implication on the analysis of cosmological observations. 

3. Distance in the inhomogeneous universe 

The Sachs e quation for the angular diameter distance is 
(ISachs 19611) 



d^ 

ds2 



(4) 



where Da is the angular diameter distance, a is the 
shear and vanishes in the Friedmann limit, Rap is the 
Ricci tensor and Ra^k'^k^ = Tajsk^k^. In the comov- 
ing and synchronous coordinates, for pressure-less mat- 
ter Raf3k°'k'^ = K,pk^k° (k = Stt and G = c = 1). Since 
the above relation depends linearly on p one could naively 
presume that since {5p) — then the average distance 
should be the same as in the unperturbed case, i.e. as in 
the Friedmann model. 
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We consider a simple example that shows otherwise: 
we assume that a ~ and Ra^k^k^ is constanl^. For 
the homogeneous model Eq. (|4]) has the same form as the 
unperturbed harmonic oscillator in Eq. (jT]). However, if 
we add density perturbations, i.e. 

p = po + Sp 

then even if {5p) ~ say Sp ^ smujx, the corresponding 
harmonic oscillator is no longer given by Eq. ^ but by 
Eq. Hence the solution to Eq. ^ for the perturbed 
universe does not oscillate about the unperturbed solu- 
tion, even if (Sp) = 0. 

How does this relate to the real Universe? Is the effect 
large ? 

The real Universe on small scales is highly inhomo- 
geneous with cosmic structures such as galaxy clusters, 
super clusters, and voids. Co s mic voids are of size 20-50 
Mpc (|Hovle fc Voeelev 20041: I Jones et al. 2009^ . thus to 
test the effect of inhomogeneities on the distance-redshift 
relations we consider models with similar scales of in- 
homogeneity. Two types of models are considered: the 
non-linear Swiss-Cheese model and the linear perturba- 
tive framework. The density profile along the line of sight, 
for both models, is presented in Fig. [2] We consider three 
different methods for calculating the distance: 

1. Fully non- linear model: 

The distance will be calculated directly from Eq. 
The evolution of matter and null geodesies will be cal- 
culated within the Szekeres Swiss-Cheese model. The 
details are presented in the Appendix. 

2. Lensing approximation: 

Within the linear perturbative scheme the distance is 



Da{z)^Da{1 + Sd), 



(5) 



where Da is the angular-diameter distance in the ho- 
mogeneous universe (see Eq. ([S])). The most general 
form for 5d was derived by Pyne & Birkinshaw (2004), 
Bonvin, Durrer & Gasparini (2006), and Hui & Greene 
(2006). Excluding the contribution from the motion of 
the observer and source, and taking the leading term, 
5d reduces to 



to 
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Fig. 2. Density fluctuations along the line of sight within 
the lensing and Swiss Cheese models. 

The present-day density fluctuations along the line of 
sight are exactly the same as in the Szekeres Swiss- 
Cheese model mentioned above. This allows us to 
see how large higher-order corrections are. The evolu- 
tion of 5 is calculated using the linear approximations 
(Peebles 1980) 



2^5='^6p. 
a & 



(7) 



3. Standard, homogeneous, Friedmann model: 

The distance within the homogeneous Friedmann 
model is 



Da[z) 



and 



(8) 



^ nil ^^{\ + zf + r!fc(l + zf + ^DE 



X exp 



3/ i±^dz' 
1 + ^' 



(9) 



where w(z) is the dark energy equation of state, and 
is parametrized as 



dx 



Xe 



xV20(x), 



(6) 



w{z) = Wo + Wl 



1 + z 



(10) 



where x is the comoving coordinate dx = dz/H{z) 
and (f> is the gravitational potential that relates to 
the density perturbations Sp via the Poisson equation 
V^(/) = ^^a^5p. Equation (|6]) is equivalent to the con- 
vergence in the lensing approximation and is known as 
the Born approximation. 



For the homogeneous case, we have ii^afe^fc" ~ a ^, so for 
small distance it changes relatively slowly. Even at distance 
around 1 Gpc, for the standard cosmological model, a(t) — 
0.8ao (where oo is the scale factor at the current instant). 



4. Evolution of dark energy 

If dark energy is some kind of dynamical field, then its 
energy density evolves. Unfortunately, evolving dark en- 
ergy cannot be modelled using the Szekeres model. The 
Szekeres model is the solution of the Einstein equations 
with a dust source -I- the cosmological constant, and hence 
cannot be used to model systems with evolving pressure. 
However, evolving dark energy can be modelled using the 
lensing approximation, but only in the case when dark 
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Fig. 3. The Szekeres Swiss-Cheese model (solid line), the 
ACDM model (dashed line) , and future observations with 
2% precision scattered around the Swiss-Cheese model. 
Upper panel presents Da{z) and lower panel presents the 
H{z) in the Swiss-Cheese and ACDM model. Since the av- 
erage H{z) in the Swiss-Cheese model is almost the same 
as in the background model (see (|A.13p ). there is no visible 
difference between these two models. 

energy is homogeneous. This is because Eq. ([5]) assumes 
that the expansion rate, which enters via Da, behaves as 
in the homogeneous model. 

When dark energy is a dynamical fluid, then two cases 
need to be considered: a clustering dark energy and non- 
clustering dark energy. If dark energy behaves as a rel- 
ativistic fluid (like radiation) then its perturbations are 
described only by decaying modes (Mukhanov 2005). In 
this case, dark energy should remain almost homogeneous. 

However, if dark energy clusters (like matter) then its 
preset-day distribution can be as inhomogeneous as mat- 
ter distribution. Then the lensing approximation [i.e. Eq. 
([S])] cannot be applied. This is because when dark energy is 
inhomogeneous then the average of the expansion rate will 
not be the same as the expansion rate of the background 
model, i.e. (H) ^ Hq. This is a consequence of gradients 
in the pressure. If dark energy is inhomogeneous and pres- 
sure gradients are non-negligible then the goo component 



of the metric must be position dependent. We denote this 
by a, i.e. goo = a^(t, r) and for simplicity consider the 
spherical case. Then (Bolejko & Lasky 2008) 



where e — pde + Pmat, Pde is the energy density of dark en- 
ergy, and Pmat is the energy density of matter. To estimate 
a we assume Wde ~ — 1 which implies that e -I- p Pmat 
and —p' « p'j^^. We also assume that the distribution of 
dark energy closely follows the matter distribution, and 
that Pde ~ Spmat which imply that 

a fa ao { — ] « 1 + 3(5, 

\PoJ 

where the subscript mat has been suppressed and the 
subscript 'o' denotes the 'observer's position', i.e. ao = 
a{r — 0) = 1, po — Pmat{f — 0). The above approxima- 
tion shows that inhomogeneities in dark energy will have 
a non-negligible effect on the dynamics of the Universe. 
Also, since a ^ 1 the result of averaging is different than 
when a — 1 (i.e. when the pressure gradients are negligi- 
ble). For example, Eq. (|A.13|) becomes 

r r 

{H)~y J driR^RYa ~y J dr{R^RY{l + 3S). 



Since d is non-zero, a deviates from unity and therefore 

W^^ = Ho- 

This feature was also observed within the second-order 
perturbative scheme. Clarkson, Ananda & Larena (2009) 
showed that the average quantities (such as the expansion 
rate) within the second-order perturbative framework are 
not the same as in the unperturbed background, in other 
words (H) 7^ Ho- Thus, if dark energy can cluster and is 
as inhomogeneous as the present-day matter distribution, 
then not only the distance fails to oscillate about the un- 
perturbed value but also the expansion rate. Therefore, in 
the next section when studying evolving dark energy it is 
implicitly assumed that dark energy is homogeneous. 

5. Future observations and dark energy 

To test the effects of inhomogeneities on the results of 
the analysis of cosmological observations, we consider the 
following example. We assume that we have data from fu- 
ture, very precise measurements of Da and H - say ten 
measurements of the baryon acoustic oscillations (BAO) 
at redshifts 0.2, 0.4, 0.6, 1.8, 2, each measurement hav- 
ing 2% errors. We then take the Szekeres model and the 
lensing approximation, and in each case generate mock 
observations (ten BAO data points) . An example of these 
data is presented in Fig. |31 where generated observations 
are scattered around the Szekeres model. We later analyse 
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Fig. 4. 68% and 95% CL for the parameters wq and wi inferred from future (mock) observations. Solid lines present 
the constraints from H(z) data, dotted line from Da{z) data, the cross marks the true value of wo and wi used to 
generate the data. The misalignment of distance-constraints is a result of not taking into account inhomogeneities. 
Left: The results of analysis within the non-linear Swiss-Cheese model {wq = — 1 and wi =0). Middle: The results 
of analysis within the linear approximation scheme {wq = — 1 and wi =0). Right: The results of analysis within the 
linear approximation scheme (wq = —0.95 and wi — —0.5). 



these mock observations using the standard methods, i.e. 
to find the equation of state of dark energy {wq and wi) 
we analyse the observations by fitting Eqs. ([8]) and dH) 
to the data. Two scenarios are considered: non-evolving 
dark energy (i.e. the cosmological constant) and an evolv- 
ing one. 

The results are presented in Fig. S) The case of non- 
evolving dark energy is presented in the left and middle 
panels. As seen if the data is analysed only using the dis- 
tance measurements, the inferred equation of state of dark 
energy differs from w = —1. This is because the inho- 
mogeneities are not taken into account in the distance- 
redshift Eq. ([5]). On the other hand, the analysis of the 
H{z) data agrees with the true equation of state of dark 
energy. This is because, within the Swiss-Cheese model 
and within the lensing approximation, the average expan- 
sion rate coincides with the unperturbed value. However, 
as pointed out in Sec. [H when dark energy is not homo- 
geneous, pressure gradients imply that (H) ^ Hq. In this 
case, if pressure gradients are large and not taken into ac- 
count then the inferred value from H{z) data may also 
differ from the true one. 

The evolving dark energy case is presented in the right 
panel of Fig. [H Only the lensing approximation is consid- 
ered here, as the Szekeres model cannot be used to de- 
scribe the evolution of dark energy. As seen, the results 
are very similar to the case of non-evolving dark energy. 

The above examples show the importance of inhomo- 
geneities, as not taking them into account can lead to 
large systematics, particularly when measuring the dark 
energy equation of state. Evolution of dark energy is a 
subtle effect, hence even small deviations of Da{z) from 
the standard case might introduce significant changes in 
the interpretation of cosmological observations. Another 
important result is that non-linear effects are important. 
As seen by comparing the left and middle panels of Fig. 
m the non-linear corrections cannot be ignored. 



6. Conclusions 

How do inhomogeneities that are present and observed in 
the Universe influence the distance-redshifl relation? 
This question has been frequently addressed since the 
papers by Kristi an fc Sachs (19 66) and Dyer & Roeder 
(1972) (see also ( Rasanen 2009f l and references therein). 
However, it is still common among cosmologists to disre- 
gard the effect of inhomogeneities. The common reasoning 
is based on the following two premises: 1) even if density 
variations are large the fluctuations of the gravitational 
potential are small and therefore the perturbation scheme 
can be applied, and 2) since perturbations are Gaussian 
they vanish after averaging and therefore they should have 
little impact on observations. 

However, as shown in this paper the distance-redshift 
relation does depend on density fluctuations (not on 
gravitational potential), and secondly, even if pertur- 
bations vanish after averaging they do modify the 
distance-redshift relation and the final result deviates 
from the homogeneous one. Even within the pertur- 
bative scheme, it is apparent that inhomogeneities do 
affect the distance-red s hift relation dBonvin et a l . 2006 : 



Vanderveld et al. 200?! : iGurzadvan fc Kocharvan 20091 ). 



Since the evolution of dark energy is a subtle effect, in- 
homogeneities cannot be ignored. 

The main aim of this paper was therefore to show 
the importance of inhomogeneous models in the studies 
of dark energy. Ignoring inhomogeneities can have seri- 
ous consequences since even small-scale inhomogeneities 
do affect the observations, especially when it comes to the 
equation of state of dark energy. This paper uses the ex- 
plicit example to show that if cosmological observations 
are analysed within the framework of homogeneous mod- 
els then matter inhomogeneities might be mistaken for 
evolving dark energy. This is a consequence of Eq. (g]), 
which is a differential equation for the angular diameter 
distance. As in the case of non-idealistic harmonic oscilla- 
tor, the perturbed solution does not oscillate around the 
unperturbed value. As trivial as it sounds, this statement 
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must be emphasized because a large number of cosmolo- 
gists seems to believe otherwise. 

The results of this paper imply that systematics due 
to inhomogeneities can distort our analysis, and may lead 
to a conclusion that dark energy evolves when in fact it is 
constant (or vice versa). Since the robust distinction be- 
tween the evolution and non-evolution of dark energy is 
the first step to understanding its nature, a proper han- 
dling of systematics is essential. Without taking into ac- 
count all systematics, the precision cosmology will not be 
an accurate cosmology (Peebles 2010). 
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Appendix A: Szekeres Swiss-Cheese model 

The metric of the Szekeres model has the form 



where 



E 



1 - k 



X- P 



S 



v_Q 

s 



1 



(A.2) 



and the prime denotes the partial derivative with respect 
to r, R' — dR/dr. The Einstein equations for the dust are 



• , 2M 1 , 



(A.3) 
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up 



2M' - 6ME'/E 
i?2(i?' - RE'/EY 



(A.4) 



where the dot denotes the partial derivative with respect 
to t, R = dR/dt, A is the cosmological constant, and M 
is an arbitrary function of r. The shear is 



1 I R' - RR'/R 



diag(0,2,-l,-l) 



3 \R'-RE'/E 
The bang time function follows from (|A.3p 

d$ 

— = t - tB{r). 



(A.5) 



(A.6) 



2M/$ + i 



A$2 



The Weyl curvature decomposed into its electric and mag- 
netic part is 



M(3i?' - RM'/M) 
3i?3(i?' - RE'/E) 



diag(0,2,-l,-l), 



(A.7) 



The constructi on of the Szekeres Swi ss-Cheese model 
closely follows (|Boleiko fc Celerier 201^ and the reader 
is referred there for more details. The functions that de- 
fine the Szekeres model were chosen to be of the form 



M = Mo + < 



f Mi£3 
Af 2 exp 







for I < Xa, 
for Xa ^ I ^ iXa 

for 3xa ^ ^ < 4a;a, 
for ^ > 4xa, 

(A.8) 

where i = r/kpc, Afo = (47rG/3c2)pfc^^ = 

nm = 0.25, iJo = 72 km s"! Mpc~i, Xa = 10^ Afi = 

2;-3M2e-i-5, Ma = -7 X 10" kpc. 

The bang time function is = 0, which means that 
the age of the universe is the same everywhere. The func- 
tion k then follows from Eq. (IA.6I) . And the functions 5*, P, 
and Q are 



5 = (5 X 10^ 
P = 1 = xo, 
Q = 1 = 2/0- 



±0.78 



(A.9) 
(A.IO) 
(A.ll) 



where -f- is for propagation from the origin [E' / E = 
0.78/(5 X 10^ + I)], and - towards the origin [E'/E = 
-0.78/(5 X 10^ €)]. As can be seen from Eqs. (|A.4j) - 
(jA.8|) . for r > 40 Mpc the considered model becomes the 
homogeneous Friedmann model, which is in this particular 
case the ACDM model. First light propagates towards the 
center, E' > 0, and after passing through the origin, E' 
becomes negative, and so on. Because the inhomogeneous 
blocks are matched to the Friedmann model the average 
density is almost the same as in the ACDM model 



-oo — oo 



(p) = ^ J dr J dp J dq 

— OO — OO 



i?2 R' - RE'/E 

E^ yr^fc 

Mo 



47r f , R^R' If, M' 





^ (A.12) 



The first equalit y is implied by the dipole not contributing 
to the average iBoleiko 200911 . The last approximation is 
exact when curvature vanishes, although in the considered 
model k < 10^^, so it is quite an accurate approximation. 
Finally, the junction conditions (matching the Szekeres to 
Friedmann model) imply that both the total mass and 
volume of the inhomogeneous patch are the same as in 
the homogeneous model. Thus, the average density of the 
Szekeres Swiss-Cheese model is almost the same as the 
background model - the ACDM model. This is also true 
for the Hubble parameter 



— oo — oo 

1 R' + 2RR'/R - 3RE'/E 
3 R' -RE'/E ~ 

47r } , R^R' ( R' ^r\ 
dr h 2- 



3V J V R' R 

^ ' 



477 



dr{R^R)' ^^=Ho 
Rf) 



(A.13) 



The redshift is obtained for null geodesic equations. 
The axially null geodesic (propagation along the axial 
axis) is 

dt R' - RE'/E 



dr ^ \/r~fc 



and the redshift relation is 



ln(l + z)^±- j dr 
or equivalently 



1 /■ , R'- RE'/E 



(A.14) 



(A.15) 



dr _ _^ 1 ^/T^ 
dz^ I + zR'-RE'/E' 
dt _ 1 R' - RE'/E 
dz ^ 1 + z R' - RE'/E' 



(A.16) 



